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Abstract. An analysis of the toroidal modes of a rotat- 
ing fluid, by means of the differential equations of motion 
is not readily tractable. A matrix representation of the 
equations in a suitable basis, however, simplifies the prob- 
lem considerably and reveals many of its intricacies. Let f2 
be the angular velocity of the star and (£, m) be the two 
integers that specify a spherical harmonic function. One 
readily finds the followings: 

1) Because of the axial symmetry of equations of mo- 
tion, all modes, including the toroidal ones, are designated 
by a definite azimuthal number m. 

2) The analysis of equations of motion in the lowest 
order of £1 shows that Coriolis forces turn the neutral 
toroidal motions of (£, m) designation of the non rotat- 
ing fluid into a sequence of oscillatory modes with fre- 
quencies 2m£l/£(£ + 1). This much is a common knowl- 
edge. One can say more, however, a) Under the Cori- 
olis forces, the eigendisplacement vectors remain purely 
toroidal and carry the identification (£, to). They remain 
decoupled from other toroidal or poloidal motions belong- 
ing to different ts. b) The eigenfrequencies quoted above 
are still degenerate. For they carry no reference to a radial 
wave number. As a result the eigendisplacement vectors, 
as far as their radial dependencies go, remain indetermi- 
nate. 

3) The analysis of equation of motion in the next 
higher order of O reveals that the forces arising from as- 
phericity of the fluid and the square of the Coriolis terms 
(in some sense) remove the radial degeneracy. The eigen- 
frequencies now carry three identifications (s,£, m), say, 
of which s is a radial eigennumber. The eigendisplace- 
ment vectors become well determined. They still remain 
zero order and purely toroidal motions with a single (£, m) 
designation. 

4) Two toroidal modes belonging to £ and £ ± 2 get 
coupled only at fi 2 order. 



5) A toroidal and a poloidal mode belonging to £ and 
£±1, respectively, get coupled but again at fl 2 order. 

Mass and mass-current multipole moments of the 
modes that are responsible for the gravitational radia- 
tion, and bulk and shear viscosities that tend to damp 
the modes are worked out in much details. 
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1. Introduction 

Recent years has witnessed a surge of interest in the small 
oscillations of rotating fluid masses. The reason for the 
excitement is the advocation by relativists that in rapidly 
rotating neutron stars the gravitational radiation drives 
the r-modes unstable, and while spinning down the star, 
may itself be amenable to detection (see recent review by 
Andersson and Kokkotas 2000 ). Nevertheless, the oscil- 
lations of rotating objects is an old problem. In the past 
few decades it has been studied by many investigators and 
from various points of view. Complexity of the problem 
arises from the fact that a fluid can support three distinct 
types of motions, derived from, say, a scalar potential, 
from a toroidal vector potential and from a poloidal vec- 
tor potential. These are the motions associated predom- 
inantly with the familiar p-, g- and toroidal- oscillations 
of the fluid. Each of these motions, in turn, can be given 
an expansion in terms of vector spherical harmonics. The 
modes of an actual star are a mixture of the three types 
mentioned above and of the different spherical harmonic 
components. Sorting out this mixture and classifying the 
modes into well defined sequences has not been an easy 
task. Moreover, and more often than not, it has not been 
realized that g- modes of spherically asymmetric config- 
urations are not apt for perturbation analysis. For, the 
low frequency tail of their spectrum is a fragile structure. 
It is driven by minute buoyancy forces and can be com- 
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pletely wiped out by almost any perturbing agent such as 
Coriolis, asphcricity and magnetic forces. 

Here we show that a matrix representation of the 
equations of motions provides a set of algebraic equa- 
tions that are much easier to cope with than their 
differential counterparts. In section 2 we write down 
the equilibrium structure and the linearized equations 
of motion of a rotating star. In section 3 we introduce 
the matrix representation of these equations. In sections 
4 and 5 we sort out the toroidal-poloidal components 
and spherical harmonic constituents of the matrices. 
In sections 6 and 7 we give an ordering of the various 
components in powers of fl 2 and sort out the equations 
of motions in various orders of magnitude. References 
and bibliographical notes relating to mode calculations 
in rotating stars are collected in section 7.4. In section 8 
we discuss the numerical results. A rotating neutron star 
can be slowed down by gravitational radiation through 
the mass and mass-current multipole moments of the 
modes. The modes, in turn, can be damped out by bulk 
and shear viscosities present in the star. The time scales 
of relevant damping mechanisms are analyzed in section 
9. Calculations of matrix elements and presentation of 
appropriate basis sets are given in the appendices. 



2. Review of rotating fluids 

2.1. Equilibrium configuration 

Let p(r, 9), p(r, 9) and U(r, 9) be the density, the pressure 
and the gravitational potential of a star rotating with the 
constant angular frequency Q about the z-axis. Equilib- 
rium condition is 



Vp + pV[U - ~fl 2 r 2 cos 2 9] = . 

For slow rotations, one obtains 

P = Po{r) + tf[p2o(r) + 0a(r)iM«»0)], 



(1) 



(2a) 



where Pi (cos 9) is a Legendre polynomial. Similar expan- 
sions exists for p and U. For a barotropic structure, p(p), 
one will have 



Po(Po), and p 2l = ^p~P2i] i = 0, 2 . 



(2b) 



is given as early as (1933) by Chandrasekhar. Further nu- 



merical values on p 2 o an d P22 may be found in Chan- 
drasekhar and Lebovitz ( 1962 ). 



2.2. Linear perturbations 

Let a mass element of the rotating fluid at position x, 
be displaced by an amount £ s (x) exp(iuj s t), where, for the 
moment, s is the collection of all indices that specify the 
displacement in question. This may include its spherical 
harmonic specifications, its radial node number, and/or 
its poloidal and toroidal nature. The Eulerian change in 
p, and U resulting from this displacement will be 



5'p = -V-(p?), 
5 s U{x) =-G j S s p{x.') I x-x 



' d 3 x' 



(3a) 
(3b) 



On the assumption that the displacement takes place 
adiabatically and the Lagrangian change in pressure is 
— (<9p/<9p)adj0V • £, one obtains 



5 s p =-(dp/dp) ad pv-e-vp-? 

= (dp/ dp) 5 s p + [dp/ dp - {dp/dp) ad ] P v ■ r , 



(3c) 



where dp /dp is the barotropic derivative of the equilib- 
rium structure. The linearized Euler equation governing 
the evolution of £ s now becomes 

+ 2ioj s pQ xf- uj s2 p£ s = , (4) 

where the linear operator W is given by 

Wf - Wp - -V P 6 s p + pV6 s U . (4a) 
P 

3. Matrix form of equations of motion 

Equation (Q) in its integro-differential form is highly com- 
plicated. We convert it into a set of linear algebraic equa- 
tions by expanding £ s in terms of a complete set of known 
basis vectors, {C"(x)}. Thus 



(5) 



where Z rs are the expansion coefficients and, as yet, are 
unknown. Associated with this basis, we define the follow- 
ing matrices 



Furthermore, Poisson's equation will give 

V 2 (C/ , 020, U22) - 47rG(p 0) P20, P22) • (2c) 

A thorough study of the structure of rotating polytropes W 



J pC*< s d 3 x, 

■H/n I pC* ■ (nxt s )d 3 x 
-in/n. J pC* x( s A, 



(6a) 

(6b) 
(6c) 



where the integration is over the volume of the star. All 
these matrices are Hermitian. Furthermore, let Z = [Z rs ] 



Y. Sobouti & V. Rezania: The r-modes of rotating fluids 



3 



be the matrix of the expansion coefficients, and to = 
[uj s 8 rs ] be the diagonal matrix of the eigenvalues. Substi- 
tuting Eq. (||) in Eq. ([I|), multiplying the resulting equa- 
tion from left by say, and integrating over the volume 
of the star gives the (q s) element of the following matrix 
equation 



WZ + 2VlCZuj - SZuj 2 = . 



(7) 



It is important to note that all factors in Eq. (]?]) are matri- 
ces and should not be commuted with each other without 
due care. Solutions of Eq. (Q) are equivalent to those of 
Eq. (H). The eigenfrequencies to in both equations are the 
same, and once the eigenmatrix Z is known, the set of 
eigendisplacement vectors {£ s } can be constructed by Eq. 

4. Partitioning into poloidal and toroidal fields 

The basis set {C r } can be divided into two poloidal and 
toroidal subsets, | Q}. The set of the eigenmodes, [£ s ], 
of Eq. (Q) in the absence of rotation also has such exact 
partitioning. Its poloidal subset comprise the commonly 
known g- and p- modes of the fluid. The toroidal subset of 
it includes those displacements of the fluid which do not 
perturb the equilibrium state of the star. For them uj = 0. 
For the sake of mathematical completeness one might say 
that uj = is a degenerate eigenfrequency and the set of 
all toroidal motions [£j ] are its eigendisplacement vectors. 

In the presence of rotation two things happen, a) Each 
known poloidal mode of the fluid acquires a small toroidal 
component, b) The neutral toroidal displacements orga- 
nize themselves, into a new sequence of modes and the 
degeneracy of uj = gets removed. Nonetheless, one may 
still partition the eigensets [£ s ] as | £t ] remembering 
that the subsets [^] and unlike the no rotation case, 
are only predominantly poloidal and toroidal, respectively. 
In view of these considerations, Eq. (pf) partitions as 



£S /-! yrs I >r ryrs 

V> — w vv ' tp ' 

St — Sp ^pt "I" St ^tt ) 



(8a) 
(8b) 



or in its matrix form and suppressing the superscripts one 
gets 



Ztp Ztt 



(8c) 



Accordingly, all matrices in Eqs. (0) partition into four 



PP, 



the block partitioned forms of all the matrices in Eq. ^ 
and extracting the tt and pt blocks of it. Thus 

tt-block of Eq.(7) : 

WuZtt + mCuZttuH - SttZuuj? 

+W tp Z pt + 2Q,C t pZ pt uj t - StpZptLOt = , (9a) 
pi-block of Eq.® : 
W pt Z tt + 2ilC pt Z u uj t — S pt Z tt uj 2 
-'rWppZpt + 2Q,CppZptLOt — SppZptLOt = . (9b) 

5. Partitioning by spherical harmonic numbers 

For a toroidal basis vector we will adopt the following 
spherical harmonic form 



=Vx[ff ( (r)17(^)] 

~ r (U 'sinr £ ' 8d >■ 
For a poloidal vector we will take 



(10a) 



i f)Y m im 

<? = -ir e (r)Yr,X rl (r)^-, X re (r)—Yn. i L0l» 



sin i? 



Appropriate ansatz for the radial function <f>, tp, and \ are 



given in appendix B. The toroidal vector ( 10a ) is obviously 
deriv ed from a radial vector potential. The poloidal vector 
( 10b ) is actually the sum of two vectors, one derived from 
a scalar potential and the ot her de rived from a toroidal 
vector potential. See Sobouti ( |l98f| ). 

In Eqs. @ we have suppressed the harmonic index 
to from £'s. For, a slowly rotating star is axially symmet- 
ric. Vectors with different values of m are not mutually 
coupled. Vectors belonging to the same to, but different 
£'s, however, are coupled. This feature entails a further 
partitioning of the basis sets into their harmonic subsets, 
0,1,2,.-.] and $,£=l,2,->>]. Correspondingly, 
each of the matrices in Eqs. (|^) partitions into blocks, 
designated by a pair of harmonic numbers (k,£), say. For 
example, the rs element of Spl, say, will be obtained from 
Eq. (^3aj) by inserting the two vectors Qp k * and Q l in that 
equation. In the following we rewrite Eqs. (|9|) taking into 
account the new partitioning. Thus, 
tt-block: 



UJ 



-W t f'z* p '* + XlCgzgut - Sgzguf } = 0, (11a) 



w£ ztf + 2nc% zfM - s% z, 



pt, tp, and tt blocks. For example, an element in the pt-block: 

pt block is obtained by inserting the two vector C and Q 

in Eq. (|a]) and carrying out the integration. 

We are not interested in the poloidal modes of Eq. 
(pa|). They are discussed in ample details and in a much 
wider scope than that of the present work in Sobouti 
(I 1980D . Here we concentrate on the toroidal modes of Eq. 
(Sb[). The required information comes from multiplying 



J pt "tt ^t 



+w£zg + XiCgzffut - Sgzguf ] = o. (Hb) 

We again emphasis that each of the factors in Eqs. (JTTJ) 
are matrices in their own right. 
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6. Expansion order and spherical harmonic 
structure of the various matrices 

Expansions of p, p, and U in powers of fl 2 results in a 
corresponding expansion of all the matrices in Eqs. (pd|). 
Moreover, having the spherical harmonics forms of Eqs. 
(|l0|), integrations over and tp dependencies in the calcu- 
lation of matrix elements can be performed analytically. 
These two tasks are carried out in appendix A. The results 
are quoted below: 



nkt qU x 



n 2 s M 



2tt( S kt 



ki 



C« = (m/£(l+l))S% t , 

Cki /~tlk* /~ikl x 
tp - ° P t — <^ot P °ke±i 



see 



(El). 



w, 



ki 
PP 



rkl 



W 0pp Sk£ 
ki 



n 2 w k fj ke± 



Wg = tfWg\ (5, 



kf: 



2tp v 
Sk,l±2) , 



<W), see (A7-|A8|), (12) 
h 8 kl±2 ), sec(|A9| -|ATc|)(13a) 
(13b) 
(13c) 
(14a) 
(14b) 
(14c) 



+ 0(n 2 ), see(A.ll), 



see(A.12), 



see(A.13), 



see(A.14) 



A subscript, 0, 2, 4 preceding the tt, tp, pt or pp designa- 
tions of the matrices indicates the order of f2 in the matrix 
in question. We also note that as f2 — * 0, u>t does the same. 
Therefore, it must have the following form 



Substituting these order of magnitude informations in Eq. 
(lib) reveals that 



rykl 
Pt 



11 A 2pt ■ 



(16a) 



For Z!ft, we have no information so far. Therefore, we 



assume the general form 

A tt — A 0tt 11 Zj 2tt ■ 



7. Expansion of equations of motion 



(16b) 



Equations ( 12 )-(|16|) allow a consistent expansion of Eqs. 
( pd| ) at VL 2 and £! 4 orders and enable one to decipher the 
information contained in them. 

7.1. Q 2 order of the tt-block 
At fl 2 order Eq. ( 11a ) gives 



ZLy 0tt Z/ 0ti °0U Zj 0tt u, 0t ~ u ■ 



(17) 



For k = £, considering the proportionality of C^ tt and S^ tt , 
Eq. (13b), one obtains 



J ott 



= 2m/£(£ + 1)1, I = unit matrix , 
undetermined at this stage. 



(18a) 
(18b) 



For k ^ £, using Eq. ( 13b ) and (18a), one has 



2m[l/fc(fc+ 1) - l)]S^Xtt = or Z& = 0.(18c) 



(|l7|). In this equation the Coriolis forces remove the de- 
generacy of the neutral motions and create a sequence of 
modes of purely toroidal nature. The new modes have a 
definite ^-symmetry, (Eq. [l^) . They are not coupled with 
toroidal motions oik ^ £ symmetry (that is Z^ t — 0) and 
with poloidal motion (that is Z^ pt = 0, see Eq. 16a). Re- 
moval of degeneracy, however, is partial. For, of the three 
designations of a standing wave in three dimensions, only 
(I, m) designations have appeared in the expression for 
a>Q t = 2m/£(£ + 1). A third designation, indicating the ra- 
dial wave number, is as yet absent. They will appear at 
higher orders of Q. 

One simplifying feature: We note that Wq 4 of Eq. (18a) 
is a constant matrix. Therefore, it will commute with 
any matrix carrying the same tl designations, such as 
Z^, Sqx, etc. This feature was used in the derivation 
of Eq.(18c) and will be used repeatedly in what follows to 
simplify the matrix manipulations. 



7.2. f2 2 order of pt -block 

Equation ( |llb| ) at order Vt 2 along with Eqs. ( |l8| ) gives 

+ (W& + 2C&?4 t )Z& =0,k = t±l. (19a) 



(15) W$hp is associated with poloidal modes of the non rotating 



fluids. In fact if £ p 's are chosen to be the exact eigenvectors 
of the non rotating star, W pp will be a diagonal matrix 
whose elements are the square of the eigenfrequencies of 
the p- and g- modes. At least i n the case of p-modes, Wq pp 



is invertible (see Sobouti 1980 for complications in the case 
of g-modes). Thus, one obtains 



yki 
J 2pt 



= -00 [w% t + c k Mzg t ,k = e±i. (i9b) 



By Eq.(8b), Eq. (19b) expresses that a toroidal mode of 
the rotating fluid of I symmetry acquires a small poloidal 
component of I ± 1 symmetry at fl 2 order. 

The case of g- modes is different. Rotation no matter 
how small cannot be treated as a perturbation on them. 
For, they are created by minute buoyancy forces and 
the low frequency tail of their spectrum gets completely 
wiped out by any other force in the medium, here 
the Coriolis forces. This, in mathematical language 
means that Wq pp for g- modes, is not invertible and 
Eq. (19b) is not applicable to them. The way out of 
the dilemma is to consider the sum of buoyancy and 
other intruding forces as an inseparable entity, without 
dividing it to large and small components. The works 
of Provost et al. (1980) and of Sobouti (1977, 1980) are 
examples of such treatments. We leave it to the experts 
in the field to decide whether the scrutiny of g- modes 
in rotating neutron stars is a crucial or an irrelevant issue. 



Let us summarize the findings so far from a pedagogical 7 3 fi 4 order of tt-block 
point of view. At fl 2 order one solves the eigenvalue Eq. 
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7.3.1. The «-subblock 

A syst emat ic extraction of f2 4 order terms of the €£-block 
of Eq. (11a) with the help of Eqs. (12)-(tL6|) and elimination 
of Z& 



rpll ytt 
1 Att Zj Qtt 



2pt term appearing in it by Eq. ( [191 
2m 



gives 



q££ 7 U , , 



0. 



£{£ + !) 

where the fourth order T-matrix is 

Tf£ = [Wut + (2C m ~ S 2 ttLOot)Lo ot 
-(W 2tp + 2Cotp^t)W^l{W 2p t + 2C m um)] 1 



(20) 



(20a) 



Equation ( fjQ[ ) is a simple eigenvalue problem. Vanishing 
of its characteristic determinant, 



rp£ 

■La 



2m 



'+!) 



Sott UJ 2t \ — : 



(20b) 



will give the non degenerate second order eigenvalues u^f, 
s = radial wave number. Once they are known, Eq. ( pp[ ) 
itself can be solved for the eigenmatrix Z^ u . We note 
that we have solved two eigenvalue problems, Eqs. ( |l7j ) 
and ( pp| ) , to remove all degeneracies of the zero frequency 
toroidal motions of the non rotating fluid. Further exten- 
sion of the analysis of equations of motion to orders higher 
than f2 4 will result is non-homogeneous algebraic equa- 
tions whose non-homogeneous terms are given in terms of 
the matrices calculated in the previous orders. 

7.3.2. The (I, £ ± 2)-subblock 



The presence of Sk,£± 2 in Eqs. (13a) and (14c) indi- 
cates that two toroidal motions belonging to £ and I ± 2 
are mutually coupled. Likewise, the presence of 6k,e±i in 



Eqs. (13c) and ( |l4bj ) shows the coupling of toroidal and 
poloidal motion of £ and £±1 symmetry. This brings in an 
additional coupling between two toroidal motions of £ and 
£±2 symmetries through the intermediary of poloidal mo- 
tions. Therefore, the only unexplored blocks of Eq. (11a) 
are those with (I, £ ± 2) designations. As in sec tion 7.3.1 
above, we extract the ft 4 order terms of Eq. (11a), but 

this time with £, I ± 2 sup erscripts, eliminate ^ 
pearing in it by Eq. (|19b|) and arrive at 



ap- 



rpi,i±2 y e±2,e±2 



where 



[2C { 



u 7 

0tt Zi 2tt 



e±2 



Q U 7 e.,t±2 /±2l, £±2 



(21) 



T±f t ±2 = [Wat + (2C 2tt - S 2 ttUot)oJat 

~ (W2tp + (£±2)(£±2 + l) Cotp) W °™ 
4m 



{W 2pt + 



(£±2)(£±2 + l) 



prCopt)] 



t,l±2 



(21a) 



In deriving this expression in two occasions we have 
substituted for co^f 2 and shifted the scalar factor 2m/ (£± 
2)(£ ±2 + 1) across the other matrices. Returning to Eq. 



( PH ) we substitute for C^ tt from Eq. (13t) and solve for 
Z 2tt . Thus, 



z: 



,£±2 



= [(C 



£±2\ 



•or 



J^ot I Wott) 1 



e±2 7 £±2,e±2 
^ott 



(22) 



Equations ©, (|l|) and (||) are all the informations con- 
tained in the S! 4 order of the tt-block. Whether Z^f t is 
non zero or otherwise is not clear at this level. To answer 
t he q uestio n on e has to go to f2 4 and tt 6 orders of Eqs. 
(lib) and (11a), respectively. This, however, will not be 
attempted here. 



7.4- Bibliographical notes 



Papaloizou and Pringlc ( 1978[ ) have studied the low fre- 
quency g- and r- modes of Eq. (^|) in an equipotential coor- 
dinate system with their applicability to the short period 
oscillations of cat aclysmic variables in mind. 

Sobouti ( |1980| ) has studied the problem primarily with 
the goal of analyzing the perturbative effects of slow ro- 
tations on p-modes and demonstrating that rotation, no 
matter how small, cannot be treated as a perturbation on 
g-modes. He argues that the g-modes are fragile structures 
and their low frequency tail of the spectrum, below the ro- 
tation frequency of the star, will be completely wiped out 
by Coriolis and asphericity forces of the star. The criterion 
for the validity of perturbation expansion is that the per- 
turbing operator should be smaller than the initial unper- 
turbed operator everywhere in the Hilbert space spanned 
by the cigenfunctions of the unperturbed operator, (Rel- 
lich, 1969] ). This condition is not met by g- modes when 
exposed to rotation, magnetic, tidal forces, etc. For, they 
have vanishingly small eigcnfrcqucncies. 



Provost et al. (1981) present an analysis of what they 
call the "quasi toroidal modes of slowly rotating stars " . 
Their work should be noted for the consistency of math- 
ematical manipulations exercised throughout the paper. 
They noted that in neutrally convective rotating stars one 
cannot have modes with predominantly toroidal motions. 
They get mixed with the neutral convective displacements. 

Lockitch and Friedman (199£) also address the hybrid 
modes with comparable toroidal and poloidal motions. 
Their work should be noted for the emphasis put on the 
(£, m) parities of the hybrid components that get coupled 
through asphericity forces. 

Yoshida and Lee ( [2000D , study Eq. (||) for a) those 
modes that are predominantly toroidal in their nature and 
b) for those that have comparable poloidal and toroidal 
components. Their latter modes are the same as those of 
Provost et al. 



(1981) 



8. Numerical result for rotating polytropes 

The bulk properties of some observed neutron stars seem 
to approximate those of a polytrope of index 1. See Ste- 
rigioulas (1998). It has become fashionable to categories 
neutron stars as stiff or soft ones depending on whether 
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their density distributions are similar to those of poly- 
tropes of index smaller than 1.5 or larger, respectively. To 
have an example of each category, sample calculations are 
given for polytropes of indices 1 and 2. 

The structure of rotating polytropes, taken from Chan- 
drasekhar (1933), is summarized in appendix C. The 
ansatz for the scalars <p se , ijj si and x st appearing in Eqs. 
( |l0| ) are given in appendix B. The required matrix ele- 
ments are reduced in appendix A. For a given N = 1, 2, . . ., 
the N x N matrices are numerically integrated. The eigen- 



values uiQ t and u>2t are calculated from Eqs. (18a) and 



(20b). Once the eigenvalues are known the various com- 



<±2,j 



ponents of the eigenmatrices Z^f t , /j 9t1 



and Z. 



;±i,. 

2pt 



are 



), (E2|)and (19b). The results are 



calculated from Eqs. ( 
given in tables 1 to 3. 

In table 1, to show the convergence of the variational 
calculations, the eigenvalues w^t are displayed for poly- 
trope of index 2, £ = m = 2, and for N — 1,2,3,4,5. 
This table should be considered as a basis for judging the 
accuracy of the numerical values. With only five varia- 
tional parameters, the first, second and third eigenvalues 
are produced with an accuracy of few parts in 10 5 , 10 4 
and 10 2 , respectively. Likewise, the numerical values in 
the remaining tables should be trusted to the same degree 
of accuracy and for the first few modes. 

In table 2, the second order eigenvalues, co^t and co- 



efficient matrices, Z^ t , Z 2 p^' , and Z^ t A ' t are given for 
polytrope of index 1 and for £ = m = 2, N = 5. In table 
3 we give the same calculations for polytrope of index 2. 
The eigenvalues are in units of \firGp. They are in agree- 
ment with those of Lindblom et al. (1999), Yoshida and 



Lee ( pOOOl) , and Morsink ( [200 1|) . Each of these authors 
have used their own technique, and different from that of 
the present paper. 

A novel feature of the present analysis is the provision 
of much detail on eigendisplacement vectors, information 
that can be profitably used to calculate any other bulk 
or local parameter of the model. For example, for modes 
belonging to £ = m = 2 one may write 



c s,2 _ f-r,1 1 7 2,2\rs 



n a EG ,3 (42)" 



an extra node in between in every curve. The general rule 
is; number of nodes for any parameter f(r, 6, tp) belong- 
ing to the radial mode number s, is s + 1. This includes 
the ever present nodes at the center and the surface of 
the star. This feature is faithfully present in all five modes 
that can be constructed from the data of table 1 and 2, 
even though we know that the numerical values for s = 4 
and 5 are only orders of magnitude. 

9. r-mode time scales 

In this section we study the dissipative effects of viscosity 
and gravitational radiation on r-modes. Quite generally 
and regardless of whether the star rotates or not, the total 
energy of an undamped mode, £(r, t), is 



E = ^ J [pkreal ' kreal + €real ' W£real]d 3 X, 

— maximum kinetic energy of the mode 



(24) 



See appendix D for proof of Eq.(24). In the presence of vis- 
cous forces and the gravitational radiation the combined 
rate of dissipation (Cutler and Lindblom, 1987) is 



d F f 

— = - (2nSa* ab Sa ab + (5a*5a)d 3 x 
at J 



wt(w t - mO) > N t (uj t - mn) M {\8D lm \ z + \8J lm \ 2 ){2h) 



l>2 



where r\ and £ are the shear and bulk viscosities respec- 
tively, 

So- a b = i(Vaib + V(,^ - |£ a bV-f), shear strain, (26a) 



Sa = V • £ bulk strain, 

4ttG (£ + !)(£ + 2) 



N, = 



c 2 ^ 1 1)[(2£+ l)!!] 2 ' 



(26b) 
c = speed of light, 



5D lm = J S P r%* md * X , mass multipole moment, (27) 



(23) 



in which the first sum is the backbone of the mode and is 
of zero order. It is toroidal motion of £ = m = 2 symmetry. 
The second sum is the coupling of £ = 3, m = 2 poloidal 
motion with £ — m — 2 toroidal motion. It is a poloidal 
motion and is of second order. The third sum is the cou- 
pling of £ = 4, m = 2 toroidal motion with £ = m = 2 
toroidal motion and again is of second order. In figures 1 
to 4 we have plotted the radial behavior of ^ Q[' 2 (Z^ 2 ) rs , 

Y, r QHz!£) rs , and ErC A (4u) rs for polytropes of in- 
dices 1,2 and the first two modes, s = 1,2. The center 
and the surface are nodes in all curves. For s = 2 there is 



5 Jim = -VW + T) J (p£ + Spv)-Y* em d 3 X, 

current multipole moment, 



(28) 



=V x {vY t 



em) 



toroidal vector spherical harmonics, 



(28a) 



For a mode of the type of Eq.(23) the bulk viscous dis- 
sipation (bv) comes entirely from its poloidal component, 
as V • Q — 0. For m = £ this has the following rate 

(dE se /dt) bv = -4 2 t n 6 {Z^ t (bv)Z^ t r,k =£+1 (29) 
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where the elements of the new matrix, (bv), are 

(bv) rq = J Q[d{r^ rk )/dr - k(k + l)x rk ] 

[d{r^ qk )/dr - k{k + l) X qk ]dr/r 3 . (29a) 

The shear viscous dissipation has contributions from both 
toroidal and poloidal components of Eq.(23). To the low- 
est order in Q, however, the toroidal component has the 
dominant contribution. Thus, 



(dE^/dt) sv = -^ t if{Z^l(sv)Z«y s , 
where 

(su) r « = £(£ +1) J r)[{rd(j) ri /dr - 2(j) rl )(rd(j) qe /dr - 2<j) q£ ) 

+ (£ 2 + £-2)^ rl ^dr/r 2 ]. (30a) 



(30) 



believed that the star looses much of its energy and angu- 
lar momentum through the gravitational radiation in this 
stage. In a case study of Andersson and Kokkotas (2000), 
the rotational period increases from 2 ms to 19 ms in one 
year. Below T = 10 s K and ft 2 - 7rGp/100, the shear 
viscosity is the dominant factor in cooling down the star. 

Acknowledgements. We wish to thank Sharon Morsink for a 
careful reading of the manuscript and illuminating comments. 
We are pleased to acknowledge the referee for his valuable com- 
ments. 



Appendix A: The elements of S, C and W 
matrices 

In calculating the elements of various matrices the follow- 
ing parameters and integrals are encountered frequently: 



Dissipation due to the gravitational radiation, to the low- 
est order, comes from the current multipolc moment, 5 r Ji. 
This has the following vector elements 



21 



S r J e = iuj< ot n— I p r e+1 (j) rl dr 



(31) 



To obtain the dissipation time scales: a)we have calculated 
the integrals of Eqs.(24)-(31), numerically, b) For the bulk 
and shear viscosities we have adopted the values of Cutler 
and Lindblom, r/ = 347p 9 / 4 T" 2 g cm" 1 s" 1 , ( = 6.0 x 
10~ 59 p 2 cJ t ~ 2 T 6 g cm -1 s -1 , where T is the temperature. 
For a mode of radial node number s and £ = m = 2, we 
obtain 



1 



•(n,T) 



1 (10 9 K 
P~ V T 

' sv \ - L 



'bv 



n 2 



- ( -I—] , 



9 2 



TlGp 



(32) 



where p is the average density of the star. Here t^ v , f bv , 



and fg T are the shear viscous-, bulk viscous- and the grav- 
itational radiation- time scale, respectively, which are nor- 
malized for T = 10 9 K and f2 2 = nGp. The total dissipa- 
tion time scale is 



1 



r(Sl,T) 



E 



i 



(33) 



For a model of 1AM@ and R = 12.53km we have calcu- 
lated ?g V , Tfo, and fg r for s = 1, 2, 3, 4, 5 and displayed in 
the unit of time in tables 4 and 5. Our values for s = 1 
arc in agreement with those of Lindblom et al. ( 1999| ) and, 
Yoshida and Lee (2000). To the best of our knowledge, the 



values for s > 2 are new. 

In a newly born hot neutron star, T > 1Q 12 K the bulk 
viscosity has a dominant role in damping out the perturba- 
tions and cooling down the star. In colder stars, T < 10 10 
K and il 2 ~ irGp, the gravitational radiation is more im- 
portant than radiation shear and bulk counterparts. While 
driving the r-modes unstable, it spins down the star. It is 



Q i = [(£ 2 -m 2 )/(2£-l)(2£+l)] 1 / 2 , 



(A.l) 



cos -dY k * m Y im d£l = Qe+ih^+i + Qeh,i-i , 

dSl = smtidMip , (A.2) 

smW k * m dY em /dMn = £Q t+1 8 k>l+1 -{£+ l)QA,£-i , 

(A.3) 

[(m 2 / sin 2 $)Y* m Y £m + 3Y*J d$dY lm / 3d] costfdfl 
£(£ + 2)Q i+x 5 k , l+l + (£ 2 - l)QM M -i , (A.4) 

cos 2 W k * m Y em dn 

(Qe+i + Qi)$kt + Qi-iQt$k,i-2 + Qe+iQe+2b~kj+2 , 

(A.5) 

[(m 2 / sin 2 $)Y k * m Y em + 3Y k \J d$dY lm / dd]P 2 {cos$)d£l 

+ ^[(£ -2){£+ l)Qi-iQt8k*-2 + l{t + 3)Q e+1 Q e+2 5 k ,e + 2] ■ 

(A.6) 

The basic definitions for S, C, and W matrices are 
given in Eqs. (^|). Expansions of p, p, and U entering these 
defining integrals are in Eqs. (^). The zero and f2 2 orders 
of these variables are sufficient to carry out the analysis 
of section [z] up to f2 4 order consistently. 

Finally, the spherical harmonics forms of the basis 
toroidal and poloidal vectors are given in Eqs. (|l^). Angu- 
lar integrals entering the definition of any matrix element 
at any desired order are performed analytically. Integrals 
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(S&y = S ki £(£ + l) I p </> r «cf> s *dr, 



in radial directions are left for numerical calculations. 
The S-matrix: 

r R 

(A.7) 



(Slf t ) rs = 5 H £(t + l) / P2^ r ^ sl dr 
Jo 

+\{5mW + 3)Q £ 2 +1 + (t- W + \)Q\ - + 1)] 
+5 kli - 2 {t - + l)Qi-iQi + 5 k j+ 2 l{i + 'i)Qe+iQw} 

x / p 2 2ci> rk 4> sl dr. (A.8) 
Jo 



The C-matrix: 

(C^r = (m/£(e+l))(S^ t r, 



R 



(A.9) 



(C 2 «r = Sum / P2o¥^ sl dT 
Jo 

9 

+^rn[5 H {Qi + i + Qt - 1/9) + S k ,e- 2 Qe-iQe 

+6k4+2Qt+iQl+»] / P22<j> rk cb sl dr. (A.10) 
Jo 



t r-ikt \rs ( r^Hk \r 

V^Opt) — V^Otp) 



= -id k , t+1 £Q e+1 / poW* + (£ + 2)x r,c W t dr 
Jo 

f R 

+i5 k , e -i(e+l)Qi / P o[V k ~ (£-l)x rk W l dr. 
Jo 

(A.ll) 



The W-matrix: 



(W&r = 5 M { / Po'dpo/dpoS^poi^Sfpoi^dr/r 2 



+ / [(d Po /dpo),d - d Po /dp }po(di; rk /dr - X rk /r) 



x (d^ se /dr- X se /r)dr/r 2 }. 



(A.12) 



= -$im[5k,t+\Qi+\ + Sk,i-iQi] 

x / p~ 1 d Po /dp p22S r p k po(r)^ l dr/r 2 , (A.13) 
Jo 

where 

S?Po(r) = [(l/ P o)(d/dr)(rp r e ) ~ t{l + l) X H ] , 

(Wg e t ) rs = 9m 2 [S ki (Q 2 +1 + Qj) + Sw-iQt-tQt 

+5k,i+2Qi+iQi+2] / p l dp /dp p 2 22 (t) rk (j) sl dr/r 2 .{k.U) 
Jo 



Poloidal motions give rise to density perturbations and 
therefore to perturbations in self gravitation. In Eqs. 
QA.12D and ( |A.13| ) the latter is neglected (Cowling's ap- 
proximation) for simplicity. Otherwise there is no concep- 
tual difficulty in including another term in Eqs. ( |A.12| -13) 
to avoid this approximation. 

Finally, we note that the upper limit of all radial inte- 
grals here is the radius of non rotating star instead of that 
of the rotating one. We reproduce the argument of Sobouti 
( 198Cl| ) to show that, in most cases, the effect arising from 
this difference in the limits of integrations is of far higher 
order in f2 2 than to upset the analysis of this paper. Let 
Ai?(#) = R($) — R be the distance between two points 
with coordinate d and situated on the surfaces of rotat- 
ing and non rotating stars. Obviously AR is of il 2 order. 
In a typical error integral, J^ +AR f(r)dr, we expand f(r) 
about R and obtain 



R+AR 



1 



f(r)dr = f(R)AR+-f'(R)(AR) 2 



(A.15) 



In Eqs. (|A.7| -14) the integrands depend on a combination 
of the variables po, P20, P22, Po and p 2 2- For a star of effec- 
tive polytropic index n, at the surface po and po vanish as 
(R — r) n and (R — r) n+1 . The leading error terms of Eq. 
(A.15) are (Ai?) n+1 cx fi 2 (" +1 ) and (Ai?)"+ 2 cx n 2{n+2 \ 
respectively. The distorted quantities P201P22 andp20)P22 
tend to zero as (R— r) n+1 and (R — r) n . However, consid- 
ering the fact that these are second order quantities and 
should be multiplied by an extra factor of fl 2 wherever 
they appear, the leading error term in expressions involv- 
ing them are also of the order S! 2 (" +1 ) and Q 2 ( n+2 \ Thus, 
the largest error committed in replacing the volume of the 
rotationally distorted star by that of the non rotating one 
is of the order Q 2 ( n+1 \ a matter of no concern for the 
analysis of this paper if n > 1. 

Appendix B: Ansatz for the scalars <p se , xp sl and 

The vicinity of the center of a star is a uniform medium, 
in the sense that, as r tends to zero, p(r), p(r), U(r), 
etc. all tend to constant values. Any scalar function, a(r) 
say, associated with a wave in such a nondispersive uni- 
form and isotropic environment should satisfy the wave 
equation V 2 tr(r) + k 2 a(r) = 0, k =const. Furthermore, 
if this scalar is associated with the spherical harmonic I, 
ie. if it is of form cr(r)Y f m (i?, (p) and is finite at the ori- 
gin, should tend to zero as r . Therefore a should have 
an expansion of the form cr(r) = ^ 'YI^LqCLsT 2 ' 1 ■ This is 
how the solutions of Laplace's equation (k = 0), spherical 
Bessel function and many other hypergeometric functions 
behave. A spherical harmonic vector, belonging to £, 
quite generally can be written in terms of three scalars 

— V(ai Y™ ) +V x A, A = ra 2 Y; n I V < (ra 3 Y e m ) ,(B.l) 
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where Cj, i = 1,2,3 are scalars of the type described Appendix D: The energy of normal mode 

above. Therefore, the radial and non radial components 
of £ should have the form 



Let us take the real part of Eq. (4) and write it in the 
following form 



(B.2) W£ e + 2f4 le x n + pf re = 



(D.l) 



We take the scalar product of Eq. (D.l) by £ rc and inte- 
To ensure this behavior, it is sufficient that the scalars grate over the volume of the star. The Coriolis term gives 
(ff l , ip rl , and \ rl entering Eqs. ( |l0| ) to be proportional to no contribution. The first term considering the hermitian 



character of W gives 



We adopt the following ansatz for the 



,X rt ) = %K +A > = 0,1,2, 



(B.3) 



where 9{r) is the polytropic function and was found by 
trial and error that ensures faster convergence of the vari- 
ational calculations. The ansatz has the required r l be- 
havior at the center. Two remarks are in order here: 

1) That a power set {r e+2r 7 r = 0, 1, 2, . . .} is complete 



1 

2 
1 d 
Ydt 



The third term in Eq. (D.l 



6-c ■ W£ TC d 3 x. 

) gives 



(D.2) 



(D.3) 



for expanding any function of r that behaves as r near Equations (O) and (O) should add to zero, which after 
theorigin follows form a theorem of Weiresstraus (Relich a timc integration gives the constant total energy 



19691 , Dixit et al. |1979p . 



2) We have chosen the ansatz of Eq. (B.3) for their sim- 
plicity. They are not the most efficient ones for rapid con- 
vergence of variational calculations. The set of the asymp- 
totic expressions that helioseismogists use for eigendis- 
placement vectors in the sun and other stars would, per- 
haps, give a faster convergence of the numerical compu- 
tations, see Christensen-Dalsgaard (1998) and references 
therein. 



Appendix C: Review of rotating polytropes 



E - E kin + E pot = - / (p£re ' £ r , 



Cre-WCre)d 3 X. (D.4) 



Next we substitute for W£ TC from Eq. (D.l). After simple 
manipulations, we obtain 



E 







1 d 2 

- 2Ekia ~ld^ 



at 



(D.5) 



The structure of rotating polytropes is taken from a land- gincc thc time dep endence of & c is sinusoidal, upon taking 

mark paper of Chandrasekhar (|1933|). A summary of what the time average of Eq . the second integral vanishes 

is needed here with slight changes in his notation is as fol- anc j are kt a i n s 
lows 

Po =Pc6 n ,P20 = np c e n - 1 y Q ,p 22 = n Pc 8 n - 1 y 2 , (C.l) 

gn+l 



E = 2£ l ki n = (maximum), 



QED. 



(D.6) 



Po = Pc& n+ ,P20 = (n + l)j3 c 6» n * ,P22 = (n + l)p c 9 n ^ 2 , References 

(C.2) 



Andersson N., and Kokkotas K. D., 2000, to appear in Int. J. 
Mod. Phys. D 

where p c and p c are constants, 9 is the polytropic variable chandrasekhar S., 1933, MNRAS 93, 390 

and satisfies the Lane-Emden equation Chandrasekhar S., Lebovitz N. R., 1962, ApJ 136, 1082 

Cutler L and Lindblom L, 1987, ApJ 314, 234 
{l/r) 2 )(d/dr])(r] 2 d9/drj) = -6 n ,r) = [(n+l)p c /4irGp 2 c ]r ,(C3) Dixit V. V., Sarath B., Sobouti Y., 1980, A&A 89, 259 

Christensen-Dalsgaard J., 1998, Lecture notes on stellar oscil- 

and ^o, ^>2 satisfy the following 
{l/r?){d/dri){Tfd*o/dri) = -nO^ 1 ^ + 1 , 



lations, 



(C.4) 

(l/7 1 2 )(d/dr 1 )(7 1 2 d^ 2 /dT 1 ) = (-710"- 1 + 6/r; 2 )* 2 . (C.5) 

It should be noted that Chandrasekhar's rotation param- 
eter is 9? /2-KGp c . For the purpose of this paper we have Lockitch K. L., Friedman J. L., 1999, ApJ 521, 764 

Morsink S., 2001, Private communications. 



tittp:/ /www. obs.aau.dk/~jcd/oscilnotes/ 
Cough D. O., 1993, in Astrophysical fluid dynamics, Les 
Houches XLVII, eds. Zahn J. P., Zinn- Justin J., Elsevier, 
Amsterdam, 399. 
Lindblom L, Mendell G, and Owen B J, 1999, Phys. Rev. D. 
60, 064006 



integrated Eqs. (|C.3| -5) numerically 



10 



Y. Sobouti & V. Rezania: The r-modes of rotating fluids 



Provost J., Berthomieu G., Rocca A., 1981, A&A 94, 126 
Papalouizou J., Pringle J. E., 1978, MNRAS 182, 423 
Rellich F., 1969, Perturbation theory of eigenvalue problem, 

Gordon and Breach Sci. Publisher, New York 
Sobouti Y., 1977, A&A 55, 339 
Sobouti Y., 1980, A&A 89, 314 
Sobouti Y., 1981, A&A 100, 319 

Ste rigioulas N., 1998, Living R eviews in Relativity 1, 8, 

bttp:/ /www. livingreviews.org 
Tassoul M. 1980, ApJS 14, 469 
Yoshida S., Lee U., 2000, ApJ 529, 997 



Figure Captions: 

— Fig. 1: Radial behavior of the various components of 
the eigenfunction of Eq.(23) for £ = m = 2, s = 1, n — 
l-ErC 2 (Zou) rs , dashed curve; E^ 3 ^)" , dot- 
dashed curve; X)r C' 4 (^2pt) rS j dotted curve. Data for 
Z's are taken from the first column of Table 2. Nodes 
in all three components are at the center and surface. 

— Fig. 2: Same as Fig. 1 for £ = m = 2, s = 2, n = 1. 
Data for Z's are from the second column of Table 2. 
Note the extra node in all three components. 

— Fig. 3: Same as Fig. 1 for £ = m = 2, s = 1, n = 2. Data 
for Z's are from the first column of Table 3. Nodes are 
at the center and the surface. 

— Fig. 4: Same as Fig. 1 for I = m = 2, s = 2,n = 2. 
Data for Z's are from the second column of Table 3. 
Note the extra nodes in all components. 



Table Captions: 



— Table 1: Convergence of the variational calculations. 
For polytrope of index 2 and I = m = 2, the second 
order toroidal eigenvalues, u)^, are displayed for differ- 
ent number of variational parameters, N = 1,2, 3, 4, 5. 
They are in units of ^pnGp. A number a x 10 ±h is writ- 
ten as a ± b. 

— Table 2: Second order eigenvalues u>2t and coefficient 
matrices, , Z^ t , and Z^ t are given for m = 2, N = 
5 and the polytropic index 1. 

— Table 3: Same as table 2, for polytrope 2. 

— Table 4: Shear viscous-, bulk viscous-, and gravita- 
tional radiation- time scales in seconds are given for 
polytrope 1 and s — 1, 2, 3, 4, 5 from top to bottom. A 
number a x 10 ±b is written as a ± b. 

— Table 5: Same as table 4, for polytrope 2. 



Tabic 1 



.30169+0 








.19860+0 


.60797+0 






.16137+0 


.44344+0 


.1371+1 




.16136+0 


.43726+0 


.11361+1 


.30568+1 


.16132+0 


.43714+0 


.10106+1 


.24613+1 .52712+1 
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Table: 2 



w 2t .32203+0 .84087+0 .19790+1 .48110+1 .96363+1 

v 2,2 
^Ott 



7 3,2 
J 2pt 



74,2 
J 2tt 



-.58565+1 
.12510+2 
-.12234+2 
.64179+1 
-.15465+1 

-.39677+0 
-.20041+0 
.63624+0 
-.41126+0 
.92562-1 

.58000+1 
-.16462+2 
.21461+2 
-.14090+2 
.37608+1 



.38532+1 
-.11134+2 
.84713+1 
-.87049+0 
-.78018+0 

-.82952+0 
.98613+0 
.96043+0 
-.37431+0 
-.53391-1 

-.30934+1 
-.10961+2 
.15392+2 
.10291+2 
-.27390+1 



-.29986+1 
.82347+1 
-.34999+1 
-.73634+2 
.49263+1 

.96949-1 
-.28392+0 
-.27949+0 
.42798+0 
-.13166+0 

.27352+1 
-.10022+2 
.14924+2 
-.10129+2 
.26484+1 



-.27754+1 
.72343+1 
.18385+0 
-.15379+2 
.11366+2 

-.12394+0 
.65146+0 
-.14499+1 
.10202+1 
-.23806+0 

.30368+1 
-.11591+2 
.18085+2 
-.12810+2 
.33165+1 



.14743+2 
-.73216+2 
.14699+3 
-.13539+3 
.47334+2 

-.27257+0 
.10959+1 
-.17471+1 
.98804+0 
-.18318+0 

-.40984+1 
.18955+2 
-.35467+2 
.30614+2 
-.10091+2 



s = 1 



Table 3 



7 3,2 
Zj 2pt 



u 2t .16132+0 

7 2,2 
^Ott 

.30158+1 
-.36757+1 

.16613+1 
-.37105+0 

.14113+0 

-.14392+1 
.10270+1 
-.28306+0 
-.17515-1 
.20617-1 

-.30726+1 
.71922+1 
-.83777+1 
.51958+1 
-.13588+1 
s = 1 



.43714+0 .10106+1 .24613+1 .52712+1 



y4.2 
J 2tt 



.17867+1 
-.25829+1 
-.30333+1 
.52330+1 
-.18708+1 

-.28790-1 
-.14184+1 
.15325+1 
.71869+0 
.13678+0 

.16563+1 
.48149+1 
-.58904+1 
.36393+1 
.93612+0 
s = 2 



.11589+1 
-.20167+0 
-.85360+1 
.14253+2 
-.59865+1 

.29383+0 
-.11345+1 
.57987+0 
-.63417-1 
-.15707-1 

-.14505+1 
.43706+1 
-.56669+1 
.34767+1 
-.85413+0 
s = 3 



.82841+0 
.15548+1 
-.14096+2 
.24110+2 
-.13019+2 

.19977-1 

.15290+0 

-.95505+0 

.68145+0 

-.15094+0 

-.14874+1 
.46605+1 
-.63205+1 
.40067+1 
-.92167+0 



.10498+2 
-.45626+2 
.84738+2 
-.75082+2 
.25924+2 

-.24544+0 
.92290+0 
-.14964+1 
.73244+0 
-.11366+0 

-.21505+1 
.81784+1 
-.13491+2 
.10782+2 
-.34032+1 
s = 5 
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Table: 4 



1~bv 



1 


2.14+8 


1.20+11 


- 2.79+0 


2 


7.78+8 


1.12+11 


- 8.71-1 


3 


1.41+9 


1.04+11 


- 5.53-1 


4 


1.72+9 


9.72+10 


- 5.66-1 


5 


1.47+8 


9.06+10 


- 1.00+0 


Tabic 5 


s 




hv 


T gr 


1 


1.28+9 


2.07+10 


-3.60+0 


2 


6.89+8 


2.84+10 


-3.51+0 


3 


1.68+8 


3.60+10 


-7.52+0 


4 


1.74+9 


4.36+10 
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